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1 Introduction 



Suppose that {Z n ;n > 1} is a Galton- Watson process with each particle having probability p n 
of giving birth to n offspring. Let L stand for a random variable with this offspring distribution. 



OO 

Let m := YlnLi n Pn be the mean number of offspring per particle. Then Z n /m n is a non-negative 
martingale. Let W be the limit of Z n /m n as n — > oo. Kesten and Stigum proved in |10| that if 
<^ \ 1 < m < oo ( that is, in the supercritical case) then W is non-degenerate (i.e., not almost surely 



zero) if and only if 



oo 

£'(Llog + L)) = ^p n nlogn < oo, (1.1) 

71=1 



here, and in the rest of this paper, we use the notation that log + r = V log r for all r > 0. This 
result is usually referred to the Kesten-Stigum L log L theorem. 

In 1995, Lyons, Pemantle and Peres developed a martingale change of measure method in [20] 
to give a new proof for the Kesten-Stigum LlogL theorem for single type branching processes. 
Later this method was extended to prove the L log L theorem for multiple and general multiple 
type branching processes in [2], [14] and [19] . 

In this paper we will extend this method to supercritical branching Hunt processes and 
establish an LlogL criterion for branching Hunt processes. To review the known results and 
state our main result, we need to introduce the setup we are going to work with first. 
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In this paper E always stands for a locally compact separable metric space. We will use 
Ea '■= E U {A} to denote the one-point compactification of E. We will use B{E) and B{Ea) to 
denote the Borel a- fields on E and Ea, respectively. Bb(E)( respectively, B + {E)) will denote the 
set of all bounded (respectively, non-negative) Z3(.E)-measurable functions on E. All functions 
/ on E will be automatically extended to Ea by setting /(A) = 0. Let ~M. P (E) be the space 
of finite point measures on E, that is, measures of the form [i = Y^i=i where n = 0, 1, 2, . . . 
and Xi £ E, i = 1, . . . , n. (When n = 0, p, is the trivial zero measure.) For any function / on E 
and any measure fi € M. P (E) we use (/, /i) to denote the integral of / with respect to fi. 

We will always assume that Y = {Yt,H x ,C} is a Hunt process on E, where £ = inf{t > : 
Yt = A} is the lifetime of Y. Let {Pt,t > 0} be the transition semigroup of Y: 

P t f(x) = U x [f(Y t )] for /e 

Let m be a positive Radon measure on E with full support. {Pt,t > 0} can be extended to a 
strongly continuous semigroup on L 2 (E,m). Let {Pt,t > 0} be the semigroup on L 2 (E,m) such 
that 

/ f(x)P t g(x)m(dx) = / g(x)P t f(x)m(dx), f,g£L 2 (E,m). 
Je Je 

We will use A and A to denote the generators of the semigroups {Pt} and {Pt} on L 2 (E,m) 

respectively. 

Throughout this paper we assume that 

Assumption 1.1 (i) There exists a family of continuous strictly positive functions {p(t, •, -);t > 
0} on E x E such that for any (t, x) £ (0, oo) x E, we have 

p tf(x)= / p(t,x,y)f(y)m{dy), P t f{x)= I p(t,y,x)f(y)m(dy). 
Je Je 

(ii) The semigroups {Pt} and {Pt} are ultracontractive, that is, for any t > 0, there exists a 
constant c< > such that 

p(t, x, y) < ct for any (x, y) <E E x E. 
Consider a branching system determined by the following three parameters: 

(a) a Hunt process Y = {Yt,H x ,C} with state space E; 

(b) a nonnegative bounded measurable function (3 on E; 

(c) an offspring distribution {(p n ( x ))'^Lo: x € E}. 
Put 

oo 

lj){x,z) =^2 Pn (x)z n , Z>0. (1.2) 
n=0 

ip is the generating function for the number of offspring generated at point x. 
This branching system is characterized by the following properties: 
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(i) Each particle has a random birth and a random death time. 

(ii) Given that a particle is born at x € E, the conditional distribution of its path after birth 
is determined by II X . 

(iii) Given the path Y of a particle and given that the particle is alive at time t, its probability 
of dying in the interval [t, t + dt) is j3{Y t )dt + o(dt). 

(iv) When a particle dies at x £ E, it splits into n particles at x with probability p n (x). 

(v) The point A is a cemetery. When a particle reaches A, it stays at A for ever and there is 

no branching at A. 

We assume that the functions p n (x), n = 0, 1, • • ■ , and A(x) := ip'(x, 1) = Y2^Lo n Pn( x ) 
are bounded S(£')-measurable and that po(x) + pi(x) = on E. The last condition implies 
A{x) > 2 on E. The assumption po(x) = on E is essential for the probabilistic proof of this 
paper since we need the spine to be defined for all t > 0. The assumption pi(x) = on E is 
just for convenience as the case p\(x) > can be reduced to the case pi(x) = by changing the 
parameters (3 and ip of the branching Hunt process. 

For any c € Bb(E), we define 



e c (t) = exp f- J 



1 c(Y s )ds ) . 



Let Xt{B) be the number of particles located in B € B(E) at time t. Then {Xt,t > 0} is a 
Markov process in M P (E). This process is called a (Y, f3, ^-branching Hunt process. For any 
fi G M.p(E), let be the law of {X t ,t > 0} when Xq = jjl. Then we have 

P^exp(-f,X t ) = exp (log u t (-),fi) (1.3) 

where Ut(x) satisfies the equation 



u t (x) = n 2 



ep(t)exp(-f(Y t ))+ / ep(s)P(Y s )^(Y s ,u t _ s (Y s ))di 



for t > 0. (1.4) 



The formula (|1.4|) deals with a process started at time with one particle located at x, and 
it has a clear heuristic meaning: the first term in the brackets corresponds to the case when the 
particle is still alive at time t; the second term corresponds to the case when it dies before t. 
The formula (jl.4p implies that 

u t {x) = Ii x f [4>(Y S , u t _ s (Y s )) - u t -s(Y s )] (3(Y s )ds + U x exp(-/(Y t )) for t > (1.5) 
Jo 

(see [3 Section 2.3]). For any fx € M P (E), f € B^{E) and t > 0, we have 

P, [(/, X t )} = H^ [e (1 _ A)fj (t)f(Y t )] . (1.6) 
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Let {P t (1 A)l3 ,t > 0} be the Feynman-Kac semigroup denned by 

P^ A)f3 f(x) := U x [f(Y t ) e (1 _ A) p(t)f(Y t )] , f € B{E). 

Let {P t (1 ~ AV ,t > 0} be the dual semigroup of {P^~ A)P ,t > 0} on L 2 (E,m). 
Under Assumption 11.11 we can easily show that the semigroups 

{Pt AV3 } and {Pt A)P } 

are strongly continuous on L 2 (E,m). Moreover, P^ 1 A '^ admits a density p^~ A ^(t, x, y) that 
is jointly continuous in (cc, y) for each t > 0: 

P^ A)P f{x) = [ pV-Wy, x, y)f(y)m(dy), for every / € B+(E). 
JE 

The generators of {P^ 1 A '^} and {P^ 1 A ^} can be formally written as A + (A — 1)(3 and 
A + (A — l)/3 respectively. 

Let cr(A + (A - l)j3) and a(A + {A - l)/3) denote the spectrum of operator A + (A - l)/3 
and A + (A — l)/3, respectively. It follows from Jentzch's Theorem (Theorem V.6.6 on page 
333 of [23] ) and the strong continuity of {P^ 1 A ^} and {-P/ 1 A ^} that the common value 
Ai := supRe(<r(A + (A — l)/5)) = supRe(cr(A + (A — l)/3)) is an eigenvalue of multiplicity 1 for 
both A + (A — l)/3 and A + [A — l)/3, and that an eigenfunction <j> of A + (^4 — l)/3 associated 
with Ai can be chosen to be strictly positive a.e. on E and an eigenfunction <p of A + (^4 — l)/3 
associated with Ai can be chosen to be strictly positive a.e. on E. By Proposition 2.3] 
we know that <p and (p are strictly positive and continuous on E. We choose <p and (p so that 
J E (p{x)(p{x)m{dx) = 1. Then 

^(x) = e- Xxt P (:i ~ A)P (j){x), 4>{x) = e~ Xlt P^~ A) ^, x £ E. (1.7) 

Throughout this paper we assume that 

Assumption 1.2 Ai > 0. 

The above assumption is the condition for the supercriticality of the branching Hunt process. 
Indeed, if Ai < 0, it is easy to see that extinction occurs almost surely from the martingale Mt((p) 
defined below. 

Let £ t = a(Y s ; s <t). Note that 

is a martingale under II X , and so we can define a martingale change of measure by 



dlla 



4>(x) 



e- Xlt e (1 _ A)p {t). 



Then {Y, lit} is a conservative Markov process, and (pep is the unique invariant probability 
density for the semigroup P^ 1 A ^ , that is, for any / £ B + (E), 

«x#,)if-^/(,w*) = f /(*)*(*)«*M<fc). 
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Let p^(t, x, y) be the transition density of Y in E under H%. Then 



P*(t,x,y) 



-Ait 



<j)(x) 



£- A W{t,x,y) <f>(y). 



Throughout this paper, we assume the following 



Assumption 1.3 The semigroups {P t 



(l-A)/3 



} and {P i 



} are intrinsic ultracontractive , that 



is, for any t > there exists a constant c% such that 

p( l ~ A ^(t, x, y) < c t 4)(x)4){y), x,yeE. 

Remark 1.4 Here are some examples of Hunt processes satisfying Assumptions li.il and \1.3\ 

(1) Suppose E = D, a domain in W*, and m is the Lebesgue measure on D. If {Y, Tl x , x £ 
D} is a diffusion killed upon leaving D with generator 

A = -V • aV + b ■ V 

where (o,ij(x))ij is uniformly elliptic and bounded with a,ij,i,j = 1, • • • ,d, being bounded functions 
in C l (M. d ) such that all their first partial derivatives are bounded, and b{,i = 1, • • • , d, are 



bounded Borel functions on R rf . It was proven in and U0j that the semigroups {P [ t 



(1-A)f) 



and {Pr )P } 

are intrinsic ultracontractive when D is a bounded lipschitz domain. For more 
general conditions on D and the coefficients for {P^ }t>o an d {P^ 1 A '^} to be intrinsic 
ultracontractive, see fW/. 

(2) Suppose E = D, a bounded open set in M. d , and m is the Lebesgue measure on D. If 
{Y, U x , x € D} is a symmetric a— stable process killed upon exiting D, where < a < 2, then 



it follows from J2J/ and U5\l that the semigroups {P^ 1 A }t>Q and {P]: 1 ^ >P }t>o are intrinsic 
ultracontractive. 

(3) For examples in which E is unbounded, see fWjj and \17j . 

(4) For more examples of discontinuous Markov processes satisfying Assumptions f"Pl and 
\1.3\ we refer our readers to f!3f and the references therein. 
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It follows from [lit Theorem 2.8] that 



1 



< ce 



-vt 



<j)(x)4>{y) 

for some positive constants c and u, which is equivalent to 



x e E, 



sup 



4>{y)4>{y) 



i 



< ce 



ut 



;i.s) 



(1.9) 



Thus for any / € Bf(E) we have 



sup 



P e/) (t,x,y)f(y)m(dy)- / (j>(y)4>(y)f(y)m{dy) 



< ce 



-ut 



4>{y)4>(y)f{y)m(dy). 
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Consequently we have 

lim JE V v yj = 1 uniformly for / G Bt(E) and x £ E. 1.10 

! E <Kv)<Kv)m™<,dy) 

For any nonzero measure € M p (£'), we define 

M t (0) := e~ Al '%^ i > 0. 
{<P, A*) 

Lemma 1.5 For any nonzero measure fi € M p (£'), Mt((f>) is a nonnegative martingale under 
P^, and therefore there exists a limit M 00 ((^)) 6 [0, oo), P^-a.s. 

Proof. By the Markov property of X, (jl.6p and (jl.7p . 



This proves that {Mt (</>), i > 0} is a non- negative P^-martingale and so it has an almost sure 
limit M oa (4>) G [0, oo) as t -> oo. □ 

In this paper we are concerned with the following classical question: under what condition 
is the limit M 00 (</)) non-degenerate (that is, P^ l (M 00 ((j)) > 0) > 0)? In [I], Asmussen and 
Hering gave a criterion for M 00 (</!)) to be non-degenerate for a general class of branching Markov 
processes under regularity conditions. More precisely, it was proved in [1] that if the underlying 
Markov process Y is positive regular (see pQ for the precise definition), M oo (0) is non-degenerate 
if and only if 

f m{dy)4>{y)P 5y [(<£, X t ) log+(0, X t )] < oo for some t > 0. (1.11) 
JE 

This condition is not easy to verify since it involves the branching process X itself. It is more 
desirable to have a criterion in terms of the natural model parameters A, (3 and {p n (x)} of the 
branching process. Such a criterion is found in [1] and [9] for branching diffusions and it was 
proved that, in the case of branching diffusions on a bounded open set E C M rf with E being 
the union of finite number of bounded C 3 -domains, M^cf)) is non-degenerate if and only if 

4>{y)P(y)l{y)m{dy) < oo. (1.12) 

JE 
where 

oo 

l(x) = ^k<j){x)\og + {k<i)(x))p k {x), xeE. (1.13) 

k=2 
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The arguments of [T] and [H] are mainly analytic. 

The purpose of this paper is two-folds. First, we generalize the result above to general 
branching Hunt processes. Even in the case of branching diffusions, our main result is more 
general than the corresponding result in [9] in some aspect since our requirement on the regularity 
of the domain is very weak. Secondly, we give a more probabilistic proof of the result, using the 
spine decomposition and martingale change of measure. Our probabilistic proof is similar to the 
probabilistic proofs of [8], [18] and [20] . 

The main result of this paper can be stated as follows. 

Theorem 1.6 Suppose that {Xt]t > 0} is a (Y, /3,tp) -branching Hunt process and that Assump- 
tions 1.1—1.3 are satisfied. Then M 0O (<^) is non- degenerate under P^ for any nonzero measure 
[i € M p (E) if and only if 

4>(x)(3(x)l(x)m(dx) < oo, (1-14) 



E 



where I is defined in (|1.13p . 



It follows from the branching property that when [i £ ~M. p (E) is given by fi = Y27=l <W ra 
1,2,..., {xi\ i = 1, • • • , n} C E, we have 



M t {4>) = Y J z 



t=i 



where X\ is a branching Hunt process starting from 8 Xi , i = 1, . . . , n. If the conclusions hold for 
the cases that fx = S x , for any x € E, then the conclusions also hold for the general cases. So 
in the remainder of this paper, we assume that the initial measure is of the form \x = 5 x ,x € E, 
and Pg x will be denoted as P x . 

This paper is organized as follows. In the next section we will discuss the spine decomposition 
of branching Markov processes. The main result, Theorem 11.61 is proved in the last section. 

2 Spine decomposition 

The materials of this section are mainly taken from [8]. We also refer to [18] for some materials. 
The main reason we present the details here is to clarify some of the points in [8]. 
Let N = {1,2, ...}. We will use 

oo 

r := \J N" 

n=0 

(where N° = {0}) to describe the genealogical structure of our branching Hunt process. The 
length (or generation) \u\ of each u G N n is defined to be n. When n > 1 and u = (ui, . . . , u n ), 
we denote (ux, . . . , u n -\) by u — 1 and call it the parent of u. For each i e N and u = (u±, . . . , u n ), 
we write ui = (u±, . . . , u n , i) for the i-th child of u. More generally, for u = (u\, . . . , u n ),v = 
(t>i, . . . , v m ) € r, we will use uv to stand for the concatenation {u\, . . . , u n , v±, . . . , v m ) of u and 
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v. We will use the notation v < u to mean that v is an ancestor of u. The set of all ancestors of 
u is given by {v £ T : v<u} = {v£T: 3w£T \ {0} such that vw = u}. The notation v < u 
has the obvious meaning that either v < u or v = u. 

A subset r C r is called a Galton- Watson tree if a) € r; b) if u, v G T, then uv £ t 
implies it 6 r; c) for all u G r, there exists r u € N such that when j G N, uj G r if and only if 
1 < i < ?V We will denote the collection of Galton- Watson trees by T. Each u G r is called a 
node of r or an individual in r or just a particle. 

To fully describe the branching Hunt process X, we need to introduce the concept of marked 
Galton- Watson trees. We suppose that each individual u G r has a mark (Y u , a u , r u ) where: 

(i) a u is the lifetime of u, which determines the fission time or the death time of particle u as 
Cu = Ylv<u °v (C0 and the birth time of uasb u = J2v<u a v (b$ = 0); 

(ii) Y u : [b u , Cm) -> gives the location of u. Given l^_i(Cu-l-) an d K, (Y u ,u G [b u , Cu)) is 
the restriction to [b u , ( u ) of a copy of a Hunt process starting from Y u _i(£ M -i — ) at time 
6 U , i.e., a process with law Hy u _ 1 (^_ 1 _) shifted by 6 U . 

(iii) r u gives the number of the offspring born by u when it dies. It is distributed as P(Y U (£ U — )) = 
(Pfc(^i(C« — )))fceN which is as defined in Section 1. 

We will use (r, Y,a, r) (or simply (r, M)) to denote a marked Galton- Watson tree. We 
denote the set of all marked Galton- Watson trees by T = {(r, M) : r G T}. 

For any r G T, we can select a line of decent £ = {£0 = 0, £1, • • • }, where £ n+ i G r is an 
offspring of £ n G r, n = 0, 1, • • • . Such a genealogical line is called a spine. We will write (M, £) 
for a marked spine. We will write n G £ to mean that n = £j for some z > 0. We will use 

f={(r,y,ff,r J : £ C r G T} 

denote the set of marked trees with distinguished spines. L t = {u G r : b u < t < ( u } is the set 
of particles that are alive at time t. 

We will use Y = {Y t )t>o to denote the spatial path followed by a spine and n = (n t : t > 0) 
to denote the counting process of fission times along the spine. More precisely, Yt = Y u (t) and 
n t = \u\, if u G Lt n £. We use nodef((r, M, £)), or simply nodet(£), to denote the node in the 
spine that is alive at time t: 

node t (£) := node t ((r, M, £)) := u if it G (C\L t . 

It is clear that nodef (£ ) = £ nt . 

If f G £, then at the fission time £ v , it gives birth to r v offspring, one of which continues the 
spine whilst the others go off to create sub-trees which are copies of the original branching Hunt 
process. Let O v be the set of offspring of v except the one belonging to the spine, then for any 
j = 1, . . . , r v such that vj G O v , we will use (r, M) v , to denote the marked tree rooted at vj. 
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Now we introduce five nitrations on T that we shall use. Define 

F t := cr{[u, r u , a u , (Y u (s), s G [b u , ( u )) : u G r G T with ( u < t) and 

[u, (Y u (s), s G [feu.t]) : u G r G T with t G [6 U> C«)]} ; 
F t := cr(F t , (node s (£), s < t)); 
Qt := a(Y s : < s < t); 

Qt ■= <r{Gu (node s (0 ■ s <t), (( u , u < node t (£))); 

Qt ■= cr(Gt, (node s (£) : s <t), {Q u ,u < node t (£)), (r u : u < node* (£))). 

The nitrations Ft, Ft, Qt, and Qt were introduced in [8], while the filtration Qt is newly 
defined. It is obvious that Q t C Qt C Qt C F t . Set 7" = Ut>o-^*> ? = Ut>o-^~t> @ = Uj>o^*> 
Q = U f >o Qt and Q = (J t > Qt- 

For each i £ £, let P 31 be the measure on (F,F) such that the filtered probability space 
(T,F, (Ft)t>o, {P x )x£e) is the canonical model for X, the branching Hunt process in E. For 
detailed constructions of P x , we refer our readers to [3], [1] and [21]. As noted by Hardy and 
Harris [8] , it is convenient to consider P x as a measure on the enlarged space T, rather than on 
T ■ We shall use P x for the restriction of P x to J-^. 

We need to extend the probability measures P x to probability measures P x on (7~, J 7 ) so 
that the spine is a single genealogical line of descent chosen from the underlying tree. We will 
assume that at each fission time we make a uniform choice amongst the offspring to decide which 
line of descent continues the spine £. Then for u G r we have 



Prob(u G = 11 — • 



It is easy to see that 

Enf- 1 - 

To define P x we recall the following representation from [19] . 
Theorem 2.1 Every Ft-measurable function f can be written as 

where f u is Ft-measurable. 

We define the measure P x on by 

dP*(T,M,oy t =dn x (Y)dL^\n) n ^fo) n i n dP^T,^?). (2.16) 



10 



where L^ Y \n) is the law of the Poisson random measure n = {{cjj : i = 1, • • • , nt} : t > 0} with 
intensity (3(Yt)dt along the path of Y, U X (Y) is the law of the diffusion Y staring from x € E, 
and Pr v (y) = Ylk>2Pk(,y)I(r v =k) ^ s the probability that individual u, located at y £ E, has r„ 
offspring. 

It follows from Theorem 12.11 that for any bounded f £ J~t, 



? x (f\Tt)= P x (^2f u (r,M)I {u ^ } 
\ueL t 



J2fu(r,M)P*(l {uee \Ft 



uGL t 



I>M0lIr- 

u<=L t v<u v 



Then we have 



P x (f) = px \Yl ^ M ) II 7" ) ' for an y bounded / e ^ * > 0. (2.17) 



In particular, 



p-(r) = p-fEIIr) 



which implies is a probability measure. P x is an extension of P x onto (T, J 7 ) and for any 
bounded / G J( we have 

■ / 'tt 6 Lt ^<n^ 

The decomposition (|2.16p of P x suggests the following intuitive construction of the system 
under P x : 

(i) the root of r is at x, and the spine process Yt moves according to the measure Tl x ; 

(ii) given the trajectory Y. of the spine, the fission time of node v on the spine is distributed 
according to L@( Y >, where is the law of the Poisson random measure with intensity 

(iii) at the fission time of node v on the spine, the single spine particle is replaced by a ran- 
dom number r v of offspring with r v being distributed according to the law P(Yf v ) = 

(Pk(XCv))k>l', 

(vi) the spine is chosen uniformly from the r v offspring of v at the fission time of v; 

(v) each of the remaining r v — l particles vj £ O v gives rise to the independent subtrees (r, M)J, 
which evolve as independent subtrees determined by the probability measure P Y <v shifted 
to the time of creation. 
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Definition 2.2 Suppose that (£l,H,P) is a probability space, {Ht,t > 0} is a filtration on 
(£l,H) and that K is a sub-a -field of H. A real-valued process Ut on (f^J 7 , P) is called a 
P(-\ K) -martingale with respect to {Ht,t > 0} if (i) it is adapted to {HtVK,t > 0}; (ii) for any 
t > 0, E(\Ut\) < oo and (Hi) for any t > s, 



We also say that Ut is a martingale with respect to {rLt,t > 0}, given K. 

Lemma 2.3 Suppose that (^l,T-L,P) is a probability space, {7it,t > 0} is a filtration on (p,,H) 
and that /Ci,/C2 are two sub-a -fields of % such that K,\ C K.2- Assume that U} is a P(-| /Ci)- 
martingale with respect to {7it,t > 0}, is a P{-\K,2)-martingale with respect to {Ht,t > 0}. 
IfUl € K.2, U? G U t , and E (\U?U?\) < oo for any t > 0, then the product U}U? is a P(-\ /Ci)- 
martingale with respect to {Ut,t > 0}. 

Proof. Suppose that t > s > 0. The assumption that U} € K2 implies that G % s V K.2- 
Then 



Lemma 2.4 Suppose that n = {{Q : i = 1, ■ ■ ■ ,nt} : t > 0} is a Poisson random measure with 
intensity (3(Y t )dt along the path ofY. Then 



E{U t \H s yK) = U s , a.s. 



P{UlU?\U s WKi)= P [P(UtU?\H s V K 2 )\H S V K x ] 



= P [UlP{U?\U s V K 2 )\Us V K x ] 
= P[UlU 2 s \U s VKi] 



= U 2 s P[Ul\U s VK x ) 

= ulul 

where in the last second equality we used the assumption that U% G % s - 



□ 




is an L^ Y > -martingale with respect to the natural filtration {Ct} of u. 



Proof. First note that 




(2.19) 



which implies that L^^- Y \rj\ ) = 1. It is easy to check that 97^ 
using the Markov property of n. We omit the details. 



(i) 



is a martingale under by 



□ 
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It follows from the lemma above that we can define a measure 

L^OO by 



dm Y ) 



C t i<iH 



V Jo 



{(A-l)P)(Y s )ds 



Lemma 2.5 For any x 6 E and t > 0, we have 



pa 



TT —- — 



1. 



(2.20) 



Proof. It follows from (|2.16|) that, given Q, for each v < £ nt , 

P a: (r(? c j|g) = ^(y c j. 

Since, given (/, {r„,i; < £ nt } are independent, we have 



1. 



□ 



The following lemma corresponds to Theorems 5.4 and 5.5 in [5] which were not proved 
there. Our results are somewhat different from those stated in Theorems 5.4 and 5.5 in [HJ. 

Lemma 2.6 (1) The process 

t ] ■= II • exp (" j\(A - l)P){Y s )ds 



is a P x (-\ Q) -martingale with respect to {Tt,t > 0}. 
(2) The process 

is a P x {-\Q)- martingale with respect to {J~t,t > 0}. 



Proof. (1) For s, i > 0, by the Markov property, we have 

-t+s 



px 



'lt+s 



p 



J] A(y c J-exp 



((A-l)p)(Y r )dr 
A(? f J • exp (- J\(A - l)P)(Y r )dr^) ■ 

C \{A-l)p)(Y r+t )dr 



V<£,n t 

px 



n ^)--p(-jf 



~(1) 

»7t ex P 



((A-l)/3)(r r+t )dr P* 



n ^ 

Sn t <"«<5n t + 3 
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For fixed t > 0, given the path of Y, the collection of fission times {{( v : £ nt < v < Cn t+3 } '■ s > 0} 
is a Poisson random measure with intensity f3(Y t+s )ds, and has law LP( Yt +-\ it follows from (|2.19j) 
that 



Then we get 



n we) 

£n t <V<£n t + s 

px 



exp / ((A-l)f3)(Y r+t )dr . 



55<i) 



(i) 



(2) For s, t > 0, by the Markov property, we have 



px 



(2) 



pa 



n 



r 



5. a ^ 

(2) 



. p x 



n 



= % 

where in the last equality we used (|2.20p . Then we have 



~{2) 



~{2) 

»7t ■ 



□ 



The effect of a change of measure using the martingale rf^ will increase the fission rate 
along the spine from (3(Yt) to (A(3)(Yt). The effect of a change of measure using the martingale 
rfj. ' will change offspring distribution from P{Yq) = {pk^Xc,i))k>i to the size-biased distribution 
P( y q) = (Pk(Yti))k>l, where p k (y) is defined by 

* f \ k Pk(y) , . 1 ^ „ 



My) 



Define 



»7t W := ~TTY ex P 



(A!-(A-l)/3)(F s )d s 



t/^ (0) is a P^-martingale with respect to \Qu t > 0}, and it is also a P x -martingale with respect 
to {Ft,t > 0}, since rfj (4>) can be expressed as 

nt 3) (0) = E <^rW«W)exp f- / (A a - (A - l)P)(Y 8 )ds) I {ue0 . (2.21) 
ueL t Wo 7 



And then we define 



Vt(<t>) = IJ r v exp y~ J Q 



((A-l)f3)(Y s )d S ) xr7 t (3) (0) 



n^n 

x ^ (2) x ^ (3) (</ ) )- 



(04-l)/3)(F s )d S ) xrf 5 ^) 
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The following result corresponds to Definition 5.6 in |S]. 
Lemma 2.7 rjt((f>) is a P x -martingale with respect to Tf 

Proof. 77* is a P x {-\ £/)-martingale with respect to {J-t,t > 0}, and is a P X {-\Q)- 
martingale with respect to {Tt,t > 0}. Note that Q C Q, and rfe € G,rfp 6 .Ft for any 
t > 0. Using Lemma 12.31 rfy rff^ is a P x (-|C/)-martingale with respect to {J~t,t > 0}. Note 
that r^ 3 '((j)) € Q and rfj^rf^ € J 7 * for any i > 0. Using Lemma [2731 again, we see that rjt(4>) = 
rf t ; 77I ;Vi ) i s a -P x -martingale with respect to {Ft,i > 0}. □ 

Lemma 2.8 Mt{4>) is the projection ofrjt((j>) onto Tt, i-e., 

M t {<j>) = P*{U<t>)\?t). 

Proof. By ([2T211 . we have 

W) = E n r « e ~ Al V(*rV(r«(*))J(«e0- 

Then 

P x (?] t (<P)\Ft) = Yl e" Al V(^)"V(^(*)) II r - ^fefll^) 
= J] e- Al V(»)~V(y«(t)) = W), 

where in the second equality we used the fact that P x {I( u &L t ni)\-^t) = I fuel*) x Iltxji r v ■ D 
Now we define a probability measure Q x on (T, F) by 



which says that on Ft, 
dQ x = Vt(0)d? x 



y = Vt(<j>), 



tAr 



D 



exp {- (Ai - (A - l)f3)(Y s )dsj dU x (Y) 



<j>{x) 

xex P (- [\(A-l)P)(Y s )ds)dL^ J! PrM*) II ^\((r,M)?) 

V ° I V<t nt j: vjeO v 

dUt(Y)dL A ^(n) [J II dP *-%((^ M )l) 

dn*(y)dL^>(n) J] p rv (Y Cv ) J] 1 n dP t !%((r,M)*). 
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Thus the change of measure from P x to Q x has three effects: the spine will be changed to a 
Hunt process under Tit, its fission times will be increased and the distribution of its family sizes 
will be sized-biased. More precisely, under Q x , 

(i) the spine process moves according to the measure 11^; 

(ii) the fission times along the spine occur at an accelerated intensity (Af3)(Yt)dt; 

(iii) at the fission time of node v on the spine, the single spine particle is replaced by a 

random number r v of offspring with size-biased offspring distribution P(Y^ v ) := (pkO^Cv)) , 

V / k>l 

where p k (y) is defined by p k (y) := k = 1, 2. • • • , y € E; 

(iv) the spine is chosen uniformly from the r v particles at the fission point v; 

(v) each of the remaining r v — 1 particles vj € O v gives rise to independent subtrees (r, M) v - 
which evolve as independent subtrees determined by the probability measure P Y ^ shifted to the 
time of creation. 

We define a measure Q x on (7~, T) by 

Q x -=Q x \t- 

It follows from Theorem 6.4 in [8] and its proof that Q x is a martingale change of measure by 
the martingale Mt((j>)- 

dQ x 



dP 1 



= M t {4>). 



Theorem 2.9 (Spine decomposition) We have the following spine decomposition for the 
martingale Mt((j>): 

Q x [<P(x)M t (^)\g] = 4>(Y t )e- x ^ + £ (r u - l)0(f f J e - A ^«. (2.22) 

Proof. We first decompose the martingale 4>(x)Mt(4>) as 

<Kx)M t (<j>) = e- Al V(?t) + e~ Al * <K y «(*))- 

The individuals {u € Lt,u ^ £ nt } can be partitioned into subtrees created from fissions along 
the spines. That is, each node u < £ nt in the spine £ has given birth at time ( u to r u offspring 
among which one has been chosen as a node of the spine whilst the other r u — 1 individuals go 
off to make the subtree (r,M)J. Put 

x i= E ^w(-), t>Cu. 

v€L t ,v€(r,M)j 

(X(,t > C u ) is a (Y, f3, ?/>)-branching Hunt process with birth time C, u and staring point Y^ u . 
Then 

<f>(x)M t (4>) = e- Al V(?t) + J] ^ M^'(^(y c Je- Al ^, (2.23) 
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where 



is, conditional on Q, a P^'-martingale on the subtree (t, M)", and therefore 

P*(M t UJ (</>)|S) = 1. 
Thus taking Q x conditional expectation of ()2.23j) gives 

Q 



Q x [4>{x)M* 
which completes the proof. 



^(y t ) e -Ai*+ (r u -mY Cu )e~ x ^, 

U<£ nt 



□ 



3 Proof of the main result 

First, we give two lemmas. The first lemma is basically [6j Theorem 4.3.3]. 

Lemma 3.1 Suppose that P and Q are two probability measures on a measurable space (Q,J- a 
with filtration (J~t)t>o> such that for some nonnegative martingale Z t , 



dQ 
dP 



= Z t . 



The limit := limsup^^ Zt therefore exists and is finite almost surely under P. Furthermore, 
for any F G J 7 ^ 

Q(F)= / z 00 dP + Q(Fn{z 00 = oo}), 
Jf 



and consequently, 



(a) P(Z 0O = 0) = 1 ^ Q(Z oc = oo) = 1 

(b) J Z^dP = j Z dP ^ Q(Zoo < oo) = 1. 



Now we are going to give a lemma which is the key to the proof of Theorem 11.61 To state 
this lemma, we need some more notation. Note that under Q x , given Q, N t := r &) : 

i = 0, 1, 2, • • • ,nt — 1} : t > 0} is a Poisson point process with instant intensity measure 
{A(3)(Yt)dtdP(Yt) at time i, where for each y € E, P(y) is the size-biased probability measure 
on N defined in Lemma 12.61 To simplify notation, Q. and rc. will be denoted as Q and r^, 
respectively. 

Recall that l(x) = X^2(^( x )) ^°& + (^(x)) Pi( x )- 
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Lemma 3.2 (1) If f E <j){y)/3(y)l(y)m(dy) < oo, then 

oo 

e- Xl ^r^(Y ) < oo, Q x - a.s. 

i=0 

^ If j E 4>{y)fi(y)l{y)m{dy) = oo, then 

limsup e~ Al ^rji^(y^ i ) = oo, Q x — a.s. 

i—^oo 

Proof. (1) For any e > 0, 

oo 
i=0 

i i 

= : I + 11. 



(3.24) 



Note that (jl.9p implies that there is a constant c > such that for any t > c and any nonnegative 
measurable function / with 

1 



OO ^ 1) 



2j E 
Then, 



4>(y)4>(y)f(y)'m(dy)< / p' t> (t,x,y)f(y)m(dy)<2 <j)(y)(p(y)f(y)m(dy), x € E. 

JE JE 



(3.25) 



Q x 
Tft 



{r 1 0(S^)>e«Ci} 



_ / 



{fc</.(r s )>e«} 



(is / p^(s,x, y)m(dy) 



ds / j/(s,x, y)m{dy) 



P(V) E fe P*(i') i W(v)>e"} 
fc=2 
oo 



fc=2 



+ / / p^(s,x, y)m{dy) 



P(y)J2kp k (y)I 



{kct>(y)>e*s} 



k=2 



(3.26) 



Applying Fubini's theorem and using the assumption that A and (3 are bounded, we get 



ds / p^(s,5c, y)m(dy) 
o Je 



/5(f) E k Pk{y)hk<t>{y)>e Ea } 
k=2 

where Ci is positive constant which only depends on c. Using (|3.25j) . we get 

oo 

P(y) E k Pk{y)hk4><,y)>e^} 



< Halloo j C ds < C x , (3.27) 

J 



/ p^(s,x, y)m(dy) 

JE 



k=2 
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< 



,. OO 

2 / m(dy)4>(y)(j)(y)/3(y)^2kp k (y) 
Je k=2 



|log+[fc0(j,)] 



ds 



£ JE 



Combining (l3T26j) . (ROTj) and (jg2Ej) , we have 



< d + - / 4>(y)P(y)l(y)m(dy). 
e Je 



_ i 

Therefore, the condition J E <j)(y) f3(y)l(y)m(dy) < co implies that 

E 7 {r^(Y Ci )>e^}<~ ; Q^-a.S 

i 

for all e > 0. Then we have 



II < oo, Q x - a.s. 



Meanwhile for e < Ai 



Q*(J) = Q a 



;,)<e £Ci } 



/'CO ^ °° 

•'0 / — o 



)<e«*} 



k=2 

oo 



f 00 _ ~ ~ Zp ~ 

< nf / dte-^- £ »P(Y t )J2kp k (Y t )I 
Jo k=2 



(3.28) 



(3.29) 



{k<t>(Y t )<e^} 



< C 2 / e- (Al - £) *cft < oo, 

JO 

where in the second to the last inequality we used the assumption that (3 and A are bounded 
and C2 is a positive constant. Then we have 



I < 00, Q x - a.s. 

Combining (pT24]l . (1339]) and (pOT)]) . we see that ^=0 e _AlCl r^(y Ci ) < 00, Q* - a.s. 

(2) It is enough to prove that for any K > 1, 

limsupe- AlCi r^(y Ci ) > K, Q x - a.s. 



(3.30) 



(3.31) 



For any fixed K > 1, define j(t, y) := f3(y) ^ k kp k (y)I 



Ut / j(t,Y t )dt 
Jo 



{fc0(j/)>A"e A i*} 



dt / m(dy)p lf '(t,x,y)j(t,y) 

JE 



jti. Since for any x G E, 
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< f dt f m{dy)p*{t,x,y)P{y)A{y) < oo, 

JO J E 

we have Jq j(t, Yt)dt < oo, lit — a.s.. Note that for any T 6 (0, oo), conditional on cr(Y), 

it{z:Cie(0,T];ri>^(r Ci )-V^} 

is a Poisson random variable with intensity 7(t, a.s. Hence, to prove (|3,3ip . we just 
need to prove 



DC 



7(t, *i)cft = oo, - a.s. (3.32) 



o 



Recall the choice of the constant c > in the statement above the inequalities (|3,25p . 
Applying Fubini's theorem and (|3.25p . we get 

rco r 

UtZ^ = dt p*(t,x,y)>y(t,y)m(dy) 



E 



2 

Exchanging the order of integration in A^, we get that 



> J dtj /(t,x,y) 7 (t,y)m(dy) 

i r°° r l 

> -/ dt <j>(y)(/>(yh(t,y)m{dy) =: (3.33) 



Aoo > / <P(y)4>(y)l 3 (y) m (dy)y2 k Pk(y)i{k>K4 ) (yy 
Je *-r 



1 , ,,,, vx l0g# 
T-log(fc^(7/)) 

Ai Ai 



> 



C3 / (/}(y)^(y)P(y)m{dy)^klog(k(f)(y))p k (y)I {k>K ^ y) -i } -C4, (3.34) 



where C3 = 1/Ai and C4 = ||/3^4|| 00 (logK+cAi)/Ai. The assumption that f E 4>(y)(3{y)l{y)m(dy) 
00 says that 

4>(y)0{y)l3(y)rn(dy)^klog(k<j)(y))p k (y)I {k>(t){y) -i } = 00. 



Since 

/ 4>(y)4>(y)P(y) m ( d y) Yl k lo g( fc( ?%))pfc(y) 

< log if I (j){y)4>(y)P(y)A(y)m{dy) < 00, 
Je 

it follows from (|3.34p that 

= 00. (3.35) 

Then by (^331) . 

nf Zoo = 00. (3.36) 
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For any finite time T > 0, put Zt = Jq l(t, Y t )dt and 

A T = [ dt [ cj>(y)0(y)j(t,y)m(dy). 
Jc Je 

Then lim T ^oo Z T = Z^, and lim^oc, A T = = TitZoa = oo. 

An argument similar to (|3.33|) yields that there exists a constant C5 which is independent of 
T and sufficiently large, such that 

-^A T < UtZ T < C 5 A T . (3.37) 
^5 

By the Paley-Zygmund inequality (see, for instance, [U Ex. 1.3.8]), 

nf ( Z T > In+Zr) > (3.38) 

Now we estimate ut(Z^). 

11*2% = lit ^ l(t,Y t )dt f T \(s,Y s )ds 
Jo Jo 

= 2Ut ( j(t,Y t )dt [ i(s,Y s )ds 
Jo Jt 

f T ~ r(t+c)AT _ *T _ rT 

= 2Ut y(t,Y t )dt j(s,Y s )ds + 2Ut l(t,Y t )dt l(s,Y s )ds 

Jo Jt Jo J(t+c)AT 

= : III + IV. 
By the Markov property of Y, we have 



IV = 2Ut [ 7 (t,Y t )dtui [ l(s,Y s -t)ds 

Jo Yt J(c+t)AT 

= 2Ut [ T 7 (t,Y t )dtui [ T 7 (u + t,Y u )du 

Jo Yt J(c+t)AT-t 

< 211* jf 7 (t,Y t )dmt^ j(u,Y u )du 



T 



m* / dt*y{t,Y t ) / du / m(dy)p*(u,Y u y) 7 ( U ,y), 



where in above inequality we used the fact that j(s,y) is decreasing in s, which is obvious by 
the definition of 7. Using (|3.25|) and (|3.37p . we see that 

iv < 4A T nfz T < 4C 5 (nf z T f. 

Meanwhile, 

"(t+c)AT 



, r-i ~ nt+c)A± 

III < 2Ut / dtj(t,Y t ) / (A(3)(Y s )d 
Jo Jt 
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< 2c\\pA\\ocnt [ l(t,Y t )dt 
Jo 

= c 6 uiz T , 



where Cq = 2c\\f3A\\ 00 . Therefore, there exists a constant C > which does not depend on T, 
such that 

n£(z T ) < c(nt(z T )f. 



Then by (JOflt we get 



nf (z T > in*z r ) > (icy 1 , 



and therefore, 



ng (z^ > ±ng z T ) > ut (z T > \ntz T \ > (Acy 1 > o. 



Since limx->oo At = oo, (|3.37p and the above inequality imply that lit (Zoo = oo) > 0. Since for 
any to > 0, J Q to -y(t, Y t )dt < ||/3A||ooio < oo, the event T := {Z^ = oo} = j J °° j(t, Y t )dt = ooj 
is invariant, that is, < ~ 1 (r) = T for all t > 0, where : i > 0} are the shift operators of 
the process Y. It follows from (|1.9p and [221 Proposition X.3.9] that the process Y is Harris- 
recurrent. Thus it follows from [221 Propositions X.3.6 and X.3.10] that lit (2 M = oo) = 1, 
which is (|3.32p . And our second conclusion follows. □ 

Proof of Theorem 11.61 The proof heavily depends on the decomposition (|2.22p . 
When J E <f>(x)f3(x)l(x)m(dx) < oo, the first conclusion of Lemma 13.21 says that 



supQ^ 



)M t {<t>)\g\ < ^r M 0(y c Je- Al ^ + II^IU < oo. 



Fatou's lemma for conditional probability implies that lim inf^oo Mt((f>) < oo, Q x -a.s. The 
Radon-Nikodym derivative tells us that Mf(0) _1 is a nonnegative supermartingale under Q x 
and therefore has a finite limit Q x -&.s. So limt^oo Mt(</>) = < oo, Q^-a.s. Lemma 13.11 
implies that in this case, 

P x [M 00 {<t>)}= lim P x [M t (cf>)] = 1. 

t— too 

When J E 4>(x)/3(x)l(x)m(dx) = oo, using the second conclusion in Lemma T3.2} we can get 
under Q x , 

lim sup 4>(x)M t (4>) > lim sup (j>{YCn )( r n t ~ l)e~ Al ^ ?li = oo. 

t— ¥OC t— ¥OC 

This yields that M 0O (^)) = oo, Q x -a.s. Using Lemma I3TT1 again, we get M oo (0) = 0, P x -a.s. The 
proof is finished. □ 
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